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The  modified  phase-space  theory  of  reaction  rates  has  been  applied  to  the 
problem  of  the  three-body  recombination  and  dissociation  of  diatomic  molecules. 
The  results  illustrate  the  important  influence  of  the  weak  attractive  minimum 
in  the  third-body  interaction  potential  and  the  effect  of  barrier  penetration 
for  recombination  at  low  temperatures.  The  system  N_  +  Ar  was  used  as  a  typical 
illustrative  example,  and  good  agreement  was  obtained  between  the  theoretical 
predictions  and  the  experimental  measurements  of  the  reaction  rate  coefficients 
over  the  temperature  range  200-12000  °K. 

The  recrossing  factor  and  the  nonequi librium  factor  were  obtained  from 
Monte  Carlo  trajectory  calculations  for  states  near  the  dissociation  limit  .- 
The  trajectories  were  sampled  within  the  reaction  zone,  with  a  weight  propor¬ 
tional  to  the  equilibrium  reaction  rate,  and  numerically  integrated  in  both 
timewise  directions  to  determine  the  complete  histories  of  the  collisions.  A 
simple,  separable  function  for  the  equilibrium  transition  rate  R(C£,  Cf)  from 
initial  energy  states  to  final  energy  states  Cf,  which  could  be  characterized 
by  three  parameters ,  was  obtained  to  fit  the  numerical  data  and  was  used  to  solve 
the  steady-state  master  equation. 

Distributions  of  the  trajectories  with  respect  to  energies  and  Impact 

parameters  are  presented,  and  several  typical  trajectories  are  shown  to  illustrate 
the  important  features  of  the  collisions.  The  contribution  to  the  overall  reac¬ 
tion  rate  from  the  "complex  mechanism  was  also  obtained. 
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MONTE  CARLO  TRAJECTORY  CALCULATIONS  OF  THE  THREE-BODY  RECOMBINATION 
AND  DISSOCIATION  OF  DIATOMIC  MOLECULES 
by 

Ven  H.  Shui 

ABSTRACT 

The  modified  phase-space  theory  of  reaction  rates  has  been  applied 
to  the  problem  of  the  three-body  recombination  and  dissociation  of 
diatomic  molecules.  The  results  Illustrate  the  important  influence  of 
the  weak  attractive  minimum  in  the  third-body  interaction  potential  and 
the  effect  of  barrier  penetration  for  recombination  at  low  temperatures. 
The  system  ^  +  Ar  was  used  as  a  typical  illustrative  example,  and  good 
agreement  was  obtained  between  the  theoretical  predictions  and  the 
experimental  measurements  of  the  reaction  rate  coefficients  over  the 
temperature  range  200-12000  °K. 

The  recrossing  factor  and  the  nonequilibrium  factor  were  obtained 
from  Monte  Carlo  trajectory  calculations  for  states  near  the  dissocia¬ 
tion  limit.  The  trajectories  were  sampled  within  the  reaction  zone, 
with  a  weight  proportional  to  the  equilibrium  reaction  rate,  and  numeri¬ 
cally  integrated  in  both  timewise  directions  to  determine  the  complete 
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histories  of  the  collisions.  A  simple,  separable  function  for  the  equi¬ 
librium  transition  rate  R(e^,  c^)  from  initial  energy  states  to  final 
energy  states  ef,  which  could  be  characterized  by  three  parameters,  was 
obtained  to  fit  the  numerical  data  and  was  used  to  solve  the  steady- 
state  master  equation. 

Distributions  of  the  trajectories  with  respect  to  energies  and  impact 
parameters  are  presented,  and  several  typical  trajectories  are  showh  to 
illustrate  the  important  features  of  the  collisions.  The  contribution 
to  the  overall  reaction  rate  from  the  "complex  mechanism"  was  also 
obtained. 

I.  INTRODUCTION 

The  most  general  theoretical  investigation  of  the  problem  of  three- 
body  recombination  and  dissociation  is  that  carried  out  by  Keck  and  his 
coworkers.  This  theory  is  based  on  a  phase  space  representation  and 
is  sufficiently  general  to  include  all  previous  classical  theories  as 
special  cases.  They  applied  the  theory  to  the  three-body  recombination 
and  dissociation  of  diatomic  molecules  with  an  inert  third-body  colli¬ 
sion  partner  (including  the  II  +  H  +  H2  recombination)  and  obtained  good 
agreement  between  theory  and  experiment  in  general.  In  addition,  their 
work  illustrated  the  importance  of  such  features  as  the  weak  attractive 
van  der  Waals  forces  between  the  recombining  atoms  and  the  inert  third- 
body  atoms  for  recombination  at  low  temperatures,  and  the  effects  of 
electronically  excited  states  and  nonequilibrium  internal  distributions 
on  the  overall  recombination  and  dissociation  rates. 

The  steady-state  reaction  rate  coefficient  (to  be  compared  with 
exnerimental  measurements)  is  obtained  as  the  triple  product 
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k  -  kB(N/N.)(k/k  )  (1.1) 

u  e 

where  k  is  the  barrier  rate/'  The  values  of  recrossing  factors  (JT/IL) 

3-5 

and  nonequilibrium  factors  (k/kg)  used  in  previous  studies  were 

2  6 

deduced  using  the  results  of  Monte  Carlo  trajectory  calculations  ’ 
and  master  equation  solutions^  in  which  the  recombining  atoms  had 
similar  masses.  Also,  the  third  body  was  considered  to  be  purely  repul¬ 
sive.  Results  of  Monte  Carlo  trajectory  calculations  of  the  dissocia- 

g 

tion  of  HCi  in  Ar  have  been  presented  recently'”  which  illustrated  the 

special  features  of  collisions  involving  highly  asvmmetric  diatomic 

molecules  and  the  importance  of  rotational  motion  in  such  collisions. 

Excellent  agreement  was  obtained  between  the  theoretical  predictions 

and  the  shock-tube  measurements  of  the  dissociation  rate  coefficients 

of  HCJl  over  the  temperature  range  2500-5000  °K. 

In  this  paper  the  results  of  Monte  Carlo  trajectory  calculations 

are  presented  for  systems  which  have  a  weak  minimum  in  the  three-body 

interaction  potential,  the  system  N2  +  Ar  being  used  as  an  illustrative 

example.  We  shall  not  repeat  details  of  the  theory  and  calculations  that 

1—8 

were  already  presented  in  previous  publications.  Instead,  in  Sec¬ 
tion  II  we  shall  summarize  the  interaction  potentials  used,  and  in  Sec¬ 
tion  III  present  the  important  features  of  the  present  calculations. 

The  results  and  their  comparison  with  experiment  will  be  presented  in 
Section  IV,  and  Section  V  will  contain  our  concluding  discussions. 

II.  INTERACTION  POTENTIALS 

The  interaction  potential  (Vq)  of  the  three-body  system  (X-Y-M) 
was  separated  into  two  parts: 


4- 


v_  *  V  4-  v 
0  VXY  V3 


(2.1) 


where  Vxy  is  the  potential  of  the  recombining  atoms  (X  -  Y) ,  and  V3  is 
the  potential  due  to  the  presence  of  the  third  body  (M).  The  Morse 
potential  function 

v  ■  D  exp[-  0(r  -  r  )]  {exp[-  g(r  -  r  )]  -  2}  (2.2) 

was  used  to  represent  Vxy.  As  was  done  before,2'3  the  dumbbell  model 


and 


V3  5  VXM  for  rXM  *  rYM 


V3  =  VYM  for  rXM  »  rYM 


(2.3) 


was  assumed  for  the  potential  V3,  and  and  were  represented  by 
the  Morse  form.  Equation  (2.1). 

The  potential  parameters  for  Vxy,  V^,  and  required  in  the 

present  calculations  were  determined  in  the  same  manner  as  in  previous 
3-5 

studies.  Table  1  summarizes  these  parameters. 

III.  MONTE  CARLO  TRAJECTORY  CALCULATIONS 

The  principles  of  the  Monte  Carlo  technique  used  here  have  been 
8 

explained  previously,  and  a  more  general  discussion  may  be  found  in  the 
9-11 

literature.  Details  of  the  selection  of  initial  conditions  and  of 
the  itegration  of  equations  of  motion  are  contained  in  Reference  12. 

In  general,  total  energy  was  conserved  to  within  0.01  kT,  and  the  trajec¬ 
tories  were  followed  in  both  time-wise  directions  from  the  initial  point 
until  the  condition 

I V3 |  <  0.01  k7  and  r3  >  0  (3.1) 


was  satisfied  (r3  is  the  distance  between  the  third  body  M  and  the  center 
of  mass  of  the  molecule  XY).  However,  for  studies  involving  complexes, 
the  alternative  cutoff  condition 
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|VXY|  <  0.01  kT  and  | |  <  0.01  kT  and  r  >  0  , 
or 

|VXY|  <  0.01  kT  and  | |  <  0.01  kT  and  r2  >  0 


(3.2) 


was  also  Introduced.  Here  r^  Is  the  distance  between  the  atom  X  and 
the  center  of  mass  of  Y  -  M,  and  r2  Is  the  distance  between  the  atom 
Y  and  the  center  of  mass  of  X  -  M. 

As  was  done  In  Reference  8,  the  recrossing  correction  factor  (N/Nq) 
was  obtained  directly  from  the  results  of  the  trajectory  calculations 
by  counting.  Also,  the  "one-way"  equilibrium  transition  rate  between 
energy  states  £  u  £^  and  e  *  E^,  per  unit  volume,  per  unit  initial  and 
final  energies,  was  fitted  by  the  empirical  forms 

R(e1,  £f  )  -  AG(ef )  exp (aef  -  8e1)  (3.3a) 

and 

R(A,  e)  -  AG (e )  exp (-  aA  -  b  e)  (3.3b) 

where  A,  a,  8,  a,  b  are  constants,  G(')  is  a  relatively  slowly  varying 
function,  A  is  the  absolute  value  of  the  energy  transfer,  e  is  the  mean 
of  the  initial  and  final  energies,  and 

a  -  a  -  b/2 ,  and  8  -  a  +  b/2  .  (3.4) 

Various  moments  of  the  energy  transfer  with  respect  to  R^,  £f) 
were  used  in  the  fitting  procedure.  In  general,  the  n>-th  moment  is 
defined  as 

Qm(e)  E  C  f(e’  Ei’  Ef)  R(Ei*  Ef>  <Ei  -  Ef)"1”1  Sm  d£.  d£f 

(3.5) 

where  £  is  the  energy  of  the  surface  from  which  samples  were  taken  (defined 
in  Reference  8),  f(£,  e^,  e^)  is  a  weighting-normalizing  function,  and 


S“  -  [(£  -  £)/|£  -  E|]m 


(3.6) 
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det ermines  the  sign  of  the  odd  moments.  The  corresponding  "absolute" 
moments  |Q|ffl(e),  are  obtained  by  setting  S  =  l,and  the  moments  used 
in  this  paper  D^e) ,  correspond  to 

f(c,  e±,  Ef)  -  G(e)/G(ef)(ftv(e)  (3.7) 

where  {Qv(e)  is  the  variational  rate  for  the  surface  S(e),  defined  in 
References  2  and  8. 

Using  Equations  (3.5),  (3,7),  and  relationships  contained  in 
Reference  8  and  evaluating  the  Integral  by  Monte  Carlo  methods,  we 
obtain 


»o 

Do<0)  -  £  tO(e)/G<e(>]  (xJ/nJ)(£1  -  Sn/H0  (3.8) 

where  ■  1  if  the  trajectory  is  reactive,  Xj  ■  0  otherwise,  n^  is 

the  number  of  crossings  of  that  trajectory  in  the  reactive  direction, 

and  Nq  Is  the  total  number  of  trajectories  sampled  on  the  surface  S(e). 

This  is  the  equation  used  in  the  actual  calculation  of  the  moments  D  (e) 

m 

13 

from  data  generated  by  trajectory  calculations.  Furthermore,  follow¬ 
ing  the  procedures  used  in  Reference  8,  we  obtain 


A/<ftv(0)  -  CaB  (3.9) 

and 

Dm(e)  -  CaB  exp(e)/fw  /"  (eA  -  exp(aef  -  BEj)  s“  de±  def 

(3.10) 


which  can  be  evaluated  exactly  to  yield 


r(m  -  D!  t(Bm  -  om)/b (oB)®"1]  C  exp[(l  -  b)e] 


Plm(e)  -  . 

I  (aB/b)fn(B/o)  C  exp[(l  -  b)e]  , 


m  >  0  ; 
m  ■  0  ; 


(3.11) 
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since  only  the  first  absolute  moment  and  the  even  moments  are  required 
in  the  present  calculations.  When  a  >>  b/2,  a  good  approximation  is 

Id!  (e)  ~  C  m!  a^  m  nxp[(l  -  b)e]  .  (3.12) 

•  'm 

The  values  of  the  parameters  a,  b,  C  were  then  derived  from  these 

moments,  and  the  steady-state  master  equation  was  solved  to  yield  the 

nonequilibrium  correction  factor  k/k  . 

e 

IV.  RESULTS 


4.1  The  Equilibrium  Transition  Rate  and  Energy  Transfer  Moments 

The  low-order  energy  transfer  moments  calculated  using  Equation  (3.8) 
are  summarized  in  Table  II.  Figure  1  shows  a  plot  of  [N(e)/NQ(e)]/[l  -  exp(-kT/D)] 
versus  e  for  the  system  N0  +  Ar.  Although  the  scatter  in  the  data  are  con- 

dm 

siderable,  the  empirical  form 

N(e)/N0(e)  -  C  exp[ (1  -  b)e]  [1  -  exp(-  kT/D)]  (4.1) 

where  D  is  the  well  depth  of  the  interaction  potential  between  the  third 
body  and  one  of  the  recombining  atoms  does  correlate  the  data  with  reasona¬ 
ble  accuracy.  Thus  the  parameters  C  and  b  were  both  determined  from  this 
figure. 


The  parameter  C  is  the  high-temperature  limit  of  the  barrier  rate 

14 

recrossing  correction  factor  N(0)/Nq(0).  Recently  Keck  has  suggested 
a  new  empirical  expression 

m..  !  -  0  7  t-J !2_  (-“l--)1/2  +  -i-  )1/2] 

N„(0)  ’  lm,  +  m„  'm,  +  m_  m,  +  m„  'm„  +  m_ 


“2  ‘“1 


2  “2 


(4.2) 


which  is  able  to  correlate  all  the  data  obtained  to  date  (including 
highly  asymmetric  diatomic  molecules).  The  possible  temperature 
dependence  of  the  parameter  b  was  also  investigated,  but  the  statistics 
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were  not  good  enough  to  give  a  definite  answer,  and  the  average  value 
of  b  was  used  in  the  following  calculations.  The  values  of  C  and  b 
are  summarized  in  Table  III. 

A  possible  correlation  of  the  parameter  b  and  the  atomic  masses  of 
the  system  is  shown  in  Figure  2  where  b  is  plotted  versus  [m^/(m^  +  m2)]^2 . 
The  straight  line 

b  -  exp{-  0.25  [m3/(m1+  m2)]1/2}  (4.3) 

represents  the  data  points  quite  well.  However,  results  for  the  systems 
^2  +  or  ^2  *  are  Per^aPs  needed  to  verify  this  correlation  since 
these  systems  have  large  values  of  the  factor  +  m2)]1/2  (4.47 

and  8.09,  respectively).  The  data  plotted  in  Figure  2  include  the  sys¬ 
tems  HCJ.  +  Ar,  N2  +  Ar,  C&2  +  Ar,  02  +  Xe,  02  +  Ar,  and  I2  +  Ar.15 

The  parameter  a  was  obtained  by  taking  the  average  of  the  four  values 
of  2|d|^/D2  in  Table  II  and  using  Che  relationship 

2 1 D 1 2^2  -  aB  (4.4) 

which  follows  from  Equation  (3.11).  The  values  of  a  for  all  the  cases 
studied  are  summarized  in  Table  III  and  plotted  versus  kT/D  in  Figure  3, 
where  D  is  the  well  depth  of  the  interaction  potential  for  the  recombin¬ 
ing  molecule.  A  weak  temperature  dependence  of  a  is  definitely  observed 
in  Figure  3,  but  the  data  are  insufficient  to  determine  any  mass  depend¬ 
ence.  The  line  of  best  least-squares  fit  (solid  straight  line  in  the 
figure)  shows  scatters  within  the  statistical  errors.  Consequently,  for 
the  data  obtained  to  date,  the  parameter  a  can  be  reasonably  well  estimated 
from  the  formula 

0.232 


a  -  3.47  (kT/D) 


(4.5) 


The  average  energy  transfers  involved  in  reactive  collisions  of 

three-body  recombination  and  dissociation  processes  are  also  of  interest 

and  can  be  evaluated  directly  from  the  moments  summarized  in  Table  II. 

The  averaged  absolute  energy  transfer  per  collision,  <A>,  and  the  root- 

2  1/2 

mean-square  energy  transfer  per  collision,  <A  >  '  ,  are  plotted  in 
Figure  4  versus  kT/D  where  D  is  the  well  depth  of  the  interaction  poten¬ 
tial  of  the  molecule.  The  best  least-squares-fit  lines  (solid  lines  in 
the  figure)  correlate  the  data  reasonably  well  and  can  be  represented 
by  the  following  equations. 

<A>  -  0.10  (kT/D)-0*30 

<A2>1/2  =  0.22  (kT/D)‘°*28  (4.6) 

Tables  IVa,  IVb,  and  IVc  show  the  distributions  of  the  trajec¬ 
tories  classified  according  to  the  nature  of  their  end  states  (see 
Reference  8).  The  cases  shown  are 

(a)  N2  +  Ar  at  11500  °K  (kT/D  -  0.1),  e  »  0  ; 

(b)  N2  +  Ar  at  100  °K  (kT/D  =  0.0009),  e  -  0  ; 

(c)  Same  as  (b),  but  with  all  trajectories  that  came  from  a 

bound  N'Ar  complex  removed. 

Although  these  distributions  have  similar  general  features  as  those 
obtained  by  Keck,8  one  major  difference  is  that  the  number  of  trajec¬ 
tories  which  cross  the  energy  surface  more  than  once  in  each  direction 
increases  as  the  temperature  is  lowered.  This  produced  a  temperature 
dependence  of  the  recrossing  factor  N/KQ  although  the  total  number  of 
f  "*■  b  trajectories  did  not  vary  significantly  with  temperature.  This 
behavior  illustrates  the  effects  of  the  attractive  minimum  in  the  third- 
body  interaction  potential  at  low  temperatures.  It  also  explains  the 
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reason  why  Keck  did  not  observe  a  similar  temperature  dependence  for 
N/Nq  (he  used  purely  repulsive  third-body  Interaction  potentials). 

In  order  to  complete  the  determination  of  the  transition  kernel, 

R(c^,  ef)  and  produce  a  map  of  Its  relative  numerical  distribution,  the 
function  G(e^)  in  Equation  (3.3)  must  first  be  found.  The  function 

G(e)  -  (1  -  (-  e/5)1/2J/[l  +  (-  e/5)1/2)  (4.7) 

2 

where  6  ■  D/kT,  was  used  for  G(c^)  In  the  present  calculations  since  Keck 
found  that  it  was  able  to  represent  the  numerical  calculations  of  <ft^(e)/i£  (0) 
with  sufficient  accuracy.  The  resulting  numerical  distributions  of 
R^ei’  are  shown  in  Figure  5  for  Nj  +  Ar  at  11500  °K.  The 

figures  in  the  boxes  represent  the  mean  value  of  R(ei»  ef)/&v( 0)  for 
the  box.  Because  of  the  sampling  technique  employed,  each  value  of  e 
used  in  the  calculations  yielded  data  in  a  quadrant  that  touches  the 
diagonal  -  e^,  at  a  value  of  e  =  (e^  +  ef)/2  -  e.  Data  in  overlapping 
quadrants  agreed  within  the  statistical  errors  and  were  averaged  to  obtain 
the  results  shown.  The  total  number  of  trajectories  sampled  was  approxi¬ 
mately  three  hundred  for  each  value  of  r(e  ■  0,  -1,  -2,  -3).  The  results 
shown  were  calculated  for  £i  >  Ef  (recombination);  however,  since  R(p-1,  Rf) 
is  symmetrical  in  and  e^,  the  values  for  may  be  obtained  by 

reflection  in  the  diagonal  ■  ef.  The  maps  show  that  R(e^,  ef)  exhi¬ 
bits  a  sharp  maximum  ridge  along  which  corresponds  to  zero  energy 

transfer;  this  dies  off  rapidly  as  le^  -  efl  increases.  Substitution  of 
Equation  (3.10)  into  Equation  (3.3)  yields 

R(e1,  £f)/<FL  (0)  -  Cag  G(ef)  exp(aef  -  ge^  .  (4.8) 

Since  the  parameters  a,  g  are  simply  related  to  the  parameters  a,  b  by 
Equation  (3.4),  Equation  (4.8)  could  be  used  to  generate  the  function 
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R(e^»  The  results  are  plotted  in  the  upper  left  corner  of 

Figure  5  as  a  contour  map.  It  can  be  seen  that  the  contour  map  repro¬ 
duces  all  the  important  features  of  the  corresponding  actual  distribu¬ 
tion,  and  the  agreement  in  general  is  very  good. 

Figure  6  shows  a  plot  of  [R  (A,  e)/<#v(0)  ]/ [R(e)/<fl(0)  ]  versus  A. 

The  factor  [<ft(e)/$(0)  ]  was  included  for  normalization  purposes  so  that 
the  data  could  be  better  represented  by  a  single  line.  In  the  plot, 
the  solid  line  passing  through  the  bulk  of  the  data  points  has  a  slope 
equal  to  the  parameter  a  as  determined  above  (see  Table  III).  It  can 
be  seen  that  for  A  <  0.5,  the  data  points  deviate  considerably  from 
this  line;  however,  the  region  around  A  *  1  has  the  most  important 
contribution  in  the  determination  of  the  nonequilibrium  correction 
factor,  and  the  line  fits  the  data  very  well  in  this  region. 

It  can  therefore  be  concluded  that  the  functional  form  for  R(e^,  e^) 
and  the  associated  parameters  determined  above  (Equation  (3.3)  and 
Table  HI)  lead  to  a  very  good  representation  of  the  transition  kernel. 

At  the  same  time,  it  is  still  simple  enough  so  that  the  steady-state 
master  equation  can  be  solved  analytically  without  the  usual  reliance 
on  a  great  deal  of  numerical  computation. 

The  nonequilibrium  correction  factors  were  then  calculated  as  was 
done  in  Reference  8* 

4.2  Recombination  and  Dissociation  Rate  Coefficients 
(1)  The  "Chaperon"  Mechanism 

One  of  the  consequences  of  including  an  attractive  well  in  the 
third-body  interaction  potential  is  that  it  provides  a  "chaperon" 
mechanism  (or  "complex"  mechanism)  in  which  long-lived  X*M  complexes 


serve  as  an  intermediate  species  in  the  reaction  process,  in  addition 
to  the  usual  "energy  transfer"  mechanism.  Tince  the  trajectory  calcula¬ 
tions  provide  complete  histories  of  the  collisions,  those  trajectories 
rrhich  contained  a  hound  complex  pair  at  the  beginning  (t  -*■  -  »)  0f  tpe 
collision  could  be  identified  by  calculating  the  relative  energies  of 
the  and  Y»X  pairs.  The  results  are  shown  in  Figure  7  as  a  plot  of 
‘C^’NC  versus  kT/Dri,  ,Jhere  is  the  ratio  of  the  effective  number 

of  reactive  trajectories  from  the  complex  mechanism  to  that  from  the 
non-complex  mechanism,  and  is  the  well  depth  of  the  third-body 
interaction  potential  (see  Table  I) .  The  data  plotted  include  results 
for  the  systems  N  (X  V)  +  Ar,  N  (A  V)  +  Ar,  and  CJL(X  V)  +  Ar  and 

L  i?  U  i  g 

can  be  represented  very  well  by  the  line 

VhiC  ■  °«»  <kT'nxM)'3/2  •  (4.3) 


The  complete  recombination  process  can  now  be  represented  by  the 
following  reactions: 


X  +  M 


X  +  X*M 


X  +  X  +  M 


♦ 

5 


X-H 


t 


X2  +  I!  , 

*  X2  +  M  . 


(4.10a) 


(4.10b) 

(4.10c) 


The  usual  rate  law  is 


and 


d(X2]/dt  -  -  d[X]/2dt 

-  k3[Xj[X*M]  +  k5[X]2  [M]  , 


V'Vc 


(4.11) 


-  k3[3M]/k5[X][K]  . 


(4.12) 
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Slnce  reactions  (4.10a)  in  general  proceed  at  much  higher  rates 
than  reactions  (4.10b)  and  (4.10c),  local  thermodynamic  equilibrium 
may  be  assumed  to  exist  between  X,  M,  and  X*M.  Thus 

[X][M]/[X-M]  «  Ke  -  k2/kx  (4.13) 

where  Kg  is  the  corresponding  equilibrium  constant. 

It  follows  that 

NC/NNC  “  k3/(Kek5)  *  (4.14) 

The  classical  phase  integrals  required  for  the  determination  of  K 

e 

were  evaluated  numerically  for  a  Morse  oscillator  using  the  parameters 
listed  in  Table  III.  The  results  show  that  within  the  temperature  range 
considered  (0.5  <  kT/D^  <  5)  , 

Ke  -  C(kT/DXI{)3/2  (4.15) 

to  a  very  good  approximation,  where  C  is  a  constant.  Hence,  from 

Equations  (4.9)  and  (4.14),  the  ratio  k^/k^  must  be  nearly  independent 

of  the  temperature  within  the  temperature  range  considered. 

The  problem  of  the  "complex"  mechanism  and  the  "energy  transfer" 

16*“22 

mechanism  has  been  studied  by  a  number  of  investigators,  ”  and 

the  equilibrium  constant  (Equation  4,13)  has  also  been  calculated  in 
2  3“  26 

other  studies.  For  the  purpose  of  comparing  with  all  the  availa¬ 

ble  experimental  measurements,  however,  the  results  are  still  very 
inconclusive. 

The  creation  of  the  "complex"  mechanism  is,  of  course,  by  no  means 
the  only  effect  of  the  attractive  minimum  in  the  third-body  interaction 
potential.  Apart  from  the  bound  complexes  (H  -  B  <  0),  which  include 
both  the  truly  bound  (H  <  0)  and  the  quantum  mechanicallj;  metastable 
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(0  <  H  <  B)  complexes,  one  can  define  the  quasidimers  (B  <  H  <  B  )  as 

m 

26 

was  done  by  Kim  and  Ross.  "Energy  transfer"  can  take  place  both 
during  the  repulsive  and  attractive  phases  of  the  collisions.  In  addi¬ 
tion,  the  attractive  minimum  i.ti  the  X  +  M  interaction  enhances  reaction 
rates  by  increasing  tht  density  of  representative  points  in  the  reac¬ 
tion  region  of  the  phase  space  (see  Reference  5  for  a  discussion  on 
this  topic  for  the  H  +  H  +  ^  recombination).  These  contributions  are 
all  included  in  the  classical  phase-space  theory  of  reaction  rates 
developed  by  Keck  and  his  coworkers. 

27 

The  recent  work  of  Clarke  and  Burns  shows  some  interesting 
features  of  the  "complex"  mechanism  of  the  recombination  of  bromine. 
However,  their  comparisons  and  conclusions  are  meaningful  only  when 
the  overall  recombination  rate  is  dominated  by  the  complex  mechanism. 
The  work  of  Keck  and  his  coworkers^"  ^  and  the  present  results  clearly 
show  that  this  condition  is  satisfied  only  when  kT  <,  D^.  At  higher 
temperatures  the  energy  transfer  mechanism  becomes  more  and  more 
important,  and  the  overall  reaction  cannot  be  represented  by  the  com¬ 
plex  mechanism  alone. 

(2)  Barrier  Penetration 

Classically,  the  barrier  surface  provides  an  unambiguous  dividing 
surface  in  phase  space  which  separates  the  reactants  (X  +  Y  +  M)  from 
the  products  (XY  +  M)  of  the  reactions  under  consideration.  Ouantum 
mechanically,  however,  this  surface  has  "holes"  due  to  barrier  penetra¬ 
tion;  i.e.,  some  reactants  can  react  and  form  products  without  passing 
through  the  barrier  surface. 
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The  effect  of  barrier  penetration  on  reaction  rates  is  twofold. 

First,  it  increases  the  barrier  rate  coefficient  k^  by  an  amount  equal 

to  the  effective  recombination  rate  due  to  barrier  penetration.  Second, 

it  may  change  the  steady-state  distribution  of  the  population  of  energy 

states  by  supplying  extra  molecules  to  states  near  the  barrier  surface, 

resulting  in  a  change  of  the  nonequilibrium  correction  factor  (k/k  ) . 

e 

These  effects  must  be  evaluated  by  solving  a  new  set  of  master  equations 
which  include  barrier  penetration  as  additional  source  terms.  This 
completely  destroys  the  simplicity  shown  in  Section  III  and  Reference  8 
where  the  steady-state  master  equation  was  solved  analytically  for  a 
separable  kernel.  Fortunately,  for  most  practical  potential  energy 
barriers  (in  particular,  the  one  used  here),  the.  transmission  coeffi¬ 
cient  (probability  of  barrier  penetration)  at  an  energy  e  decays  extremely 
rapidly  as  e  varies  from  0  to  -  B/kT.  Thus  a  new  energy  surface  (e  =  e.  ) 

D 

may  be  defined,  such  that  the  recombination  rate  due  to  barrier  penetra¬ 
tion  at  energy  levels  below  is  negligible  compared  with  the  recombina¬ 
tion  rate  due  to  collisional  deactivation  through  the  energy  surface 
e  =  e^.  Using  this  surface  as  the  new  boundary  (instead  of  e  ■  0)  and 
solving  the  steady-state  master  equation  in  a  manner  similar  to  that  of 
Reference  8,  the  factor 

t«(eb)^(0)]Rk/ke)e_e  /(k/ke)emQ]  »  exp(-  beb)  (4.16) 

b 

was  obtaired  which  must  multiply  the  factor  k/kg  in  Equation  (1.1)  to 
include  effects  of  barrier  penetration. 

In  order  to  determine  eb,  an  effective  recombination  rate  coeffi¬ 
cient  due  to  barrier  penetration  is  defined  as 
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W  [x>'2  wr1  («.i7) 

where  ft,  (e^)  is  the  total  barrier  penetration  rate  for  e  <  e,  and  Is 
K  b 

given  by 

^p^eb^  ^  “  2a(8kT/^ry)1/,2  /_gykT  Tr(e)  exp(-e)  de  .  (4.18) 

Here  a  Is  a  collision  cross-section,  and  T.(e)  is  the  transmission 
coefficient  calculated  using  the  semi-classical  equation 

x 

Tr(e)  =  expt-  (2/h)  /^{2M[V(x)/kT  -  e]}1/2  dx]  (4.19) 

where  h  is  the  Planck's  constant,  V(x)  is  the  potential  energy,  and  x1 
and  x2  are  the  two  positions  at  which  V/kT  *  e.  Equations  (4.18)  and 
(4.19)  were  then  integrated  numerically.  The  results  for  nitrogen  are 
shown  in  figure  8  where  kp  is  plotted  versus  eb,  from  which  the  new 
surface  e  -  efe  can  be  approximately  located  with  the  condition 

W  "  •  (4.20) 

For  example,  for  nitrogen  at  100  °K,  k®  (il/Ng)  *  l(f32  cm6  molecule"2 
sec  which  gives  «  -  0.35,  and  from  Equation  (4.16)  the  correspond¬ 
ing  increase  in  rate  due  to  barrier  penetration  is  a  factor  of  1.3 
approximately.  It  should  be  pointed  out  that  because  k  varies  rapidly 
with  eb,  the  determination  of  efa  is  very  insensitive  to  slight  variations 

g 

in  O,  [M],  and  kr(N/H0).  Therefore,  high  accuracy  is  not  required  in 
estimating  these  quantities. 

(3)  Corrections  Due  to  Deviations  of  the  Morse  Form  from  True  Potentials 
The  Morse  form  (Equation  (2.2)  and  Table  I)  used  to  represent  the 
interaction  potential  was  found  to  deviate  considerably,  in  the 
long-range  asymptotic  region,  from  the  more  accurate  potential  curves 
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given  by  Gilmore.  A  plot  of  Gilmore's  potentials  versus  the  separa¬ 
tion  shows  that  the  form  r  ^  is  a  much  better  representation  of  the 

o 

true  potentials  in  this  region  (r  >  1.5  A).  The  major  effect  of 

incorporating  this  modification  is  that  it  alters  somewhat  the  position 

of  the  rotational  maximum  z 2  and  thereby  changes  the  barrier  rate  coeffi- 

2 

cient  (Equation  (2,9)  of  Reference  3)  through  the  factor  z2(z2  ~  z])  * 

E  B  B 

Denoting  the  new  barrier  rete  coefficient  by  kg,  the  ratio  kg/k^  is  included 

in  Table  V.  The  corrections  are  fairly  small,  but  they  do  introduce  an 

additional  weak  negative  temperature  dependence, 

3  + 

(4)  Contribtuion  from  the  A  £  State 

u 

The  overall  dissociation  process  for  nitrogen  may  be  described 
satisfactorily  by  the  following  set  of  reactions: 


N.  (X  1I+)  +  Ar  — ^  N  +  N  +  Ar 

L  8 


N-  (X  h+)  +  H  N,  (A  3£+)  +  K 

z  8  if  z  u 

*4 


(4.21a) 

(4.21b) 


N_ (A  V")  +  Ar  — H  +  H  +  Ar  .  (4.21c) 

2  u 

It  has  been  pointed  out  previously  (see  Reference  3  and  references 

3  + 

quoted  there)  that  the  A  £y  state  makes  significant  contributions  to 

the  overall  reaction  rates.  In  order  to  properly  account  for  these 

3  + 

contributions,  the  concentration  [N2(A  E  )]  must  be  found.  A  steady- 
state  analysis  yields 

[N2(aV)]s/[N2(A  V)]e  =  U  +  (k5/k4) [Arl/[N]}_1  (4.22) 

where  the  subscript  s  denotes  steady-state,  and  e  denotes  local  equi- 

3 

librium.  Since  the  ratio  k,-/k^  was  not  known,  the  previous  calculations  ’ 
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3  + 

simply  assumed  that  the  A  2^  state  was  in  local  thermal  equilibrium 
with  the  ground  state.  Recently  Flagan  and  Appleton2^  have  obtained 
the  expression 

k4  -  5.1  x  10"3  t“2,23  cm3  sec-1  (4.23) 

for  6000°  K  <  T  <  14000  °K,  by  matching  the  computed  peak  intensity  of 

the  nitrogen  first-positive  radiation  with  Wray's  experimental  measure- 
30 

ments.  The  result  of  the  steady-state  analysis,  Equation  (4.22),  may 
therefore  be  used  in  calculating  the  rate  coefficients.  A  plot  of  the 
ratio  [^(A  ^u)  38/ 1^2 (A  ^u)le  as  a  function  of  the  temperature  shows 
that  the  assumption  of  local  thermodynamic  equilibrium  between  the 
A-state  and  the  ground  state  is  valid  for  temperatures  k  4000  °K  (wh4ch 
covers  all  recombination  rate  measurements),  but  the  contribution  from 
the  A-state  decreases  rapidly  at  higher  temperatures. 

Taking  into  account  effects  of  barrier  penetration  and  corrections 
due  to  deviations  of  the  Morse  form  from  true  interaction  potentials  for 
nitrogen,  the  results  of  the  calculated  recombination  rate  coefficients 
are  summarized  in  Table  V  and  plotted  in  Figure  9  as  a  function  of  the  tem¬ 
perature.  The  theoretical  points  ♦  [(ky  (X  hp  +  ky(A  V)J  and  •  [kr(X  VJ)] 
have  been  connected  visually  by  full  lines  to  help  in  comparing  the  results. 

The  dashed  line  represents  the  theoretical  predictions  of  k  (X  1E+)  without 

r  8 

the  corrections  due  to  barrier  penetration  and  deviations  of  V„  from 

2 

the  Morse  form.  It  is  seen  from  the  figure  that  at  low  temperatures, 
the  present  theoretical  predictions  of  the  recombination  rate  coefficient, 

+  kr^A  2u^’  are  ln  good  a8ree,nent  with  the  experimental 
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31  32 

measurements  of  Campbell  and  Thrush  (A),  and  Clyne  and  Stedman  (0); 

33 

the  more  recent  room  temperature  measurement  of  Brennen  and  Shane  (V) 

Is  somewhat  higher  but  also  appears  to  be  within  the  combined  statisti¬ 
cal  errors.  There  still  exists  a  discrepancy  of  about  a  factor  of  6 
between  Clyne  and  Stedman' s  measurement  at  90  °K  and  the  theoretical 
prediction  at  100  °K.  However,  since  Clyne  and  Stedman  used  pure  nitro¬ 
gen  in  their  experiment,  the  different  efficiencies  of  i-I^  and  Ar  as 
third  bodies  may  be  responsible  for  this  discrepancy.  In  this  regard, 
experimental  measurements  at  T  a  100  °K  using  Ar  as  the  third  body 
should  be  very  useful. 

For  the  purpose  of  effecting  a  comparison,  the  experimental  dissocia¬ 
tion  rate  measurements  have  been  converted  into  equivalent  recombination 

rate  coefficients  through  the  equilibrium  constant.  The  results  of 

34 

Appleton,  Steinberg,  and  Liquornick  (A)  are  in  excellent  agreement  with 

35  36 

the  theoretical  predictions.  The  results  of  Byron  (B)  and  Cary  (C) 
appear  to  be  somewhat  higher,  but  in  view  of  the  extent  of  experimental 
error  and  the  uncertainties  in  the  theory,  the  agreement  is  considered 
to  be  very  good. 

V.  DISCUSSION 

A  number  of  trajectories  for  the  collision  of  Kj  +  Ar  are  shown  in 
Figure  10  where  two  different  coordinate  systems  have  been  used  in  order 
to  show  certain  special  features  of  the  collisions  more  satisfactorily. 
One  coordinate  system  has  the  argon  atom  stationary,  with  the  two  nitro¬ 
gen  atoms  moving  about  and  is  particularly  useful  for  illustrating  the 
individual  collisions  between  each  nitrogen  atom  and  the  argon  atom. 
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The  other  coordinate  system  has  the  center  of  mass  of  the  nitrogen 
molecule  stationary  with  all  three  atoms  in  motion.  This  system 
offers  a  good  visualization  of  the  Internal  relative  motion  of  the 
bound  nitrogen  molecule.  A  bound  NAr  complex  is  clearlv  illustrated 
in  Figure  10(c). 

It  should  be  pointed  out  that  for  convenience  of  representation 
in  two  dimensions,  the  trajectories  shown  in  Figure  10  have  all  three 
particles  moving  in  the  same  plane.  However,  the  calculations  were, 
in  general,  three  dimensional. 

Typical  reactive  trajectories  for  Nj  +  Ar  collisions  may  be 
divided  into  the  following  types: 

(i)  Simple  f  -*•  b  Trajectories 

The  three  atoms  start  as  free  particles  and  approach  close 
enough  for  energy  and  momentum  exchanges  of  appreciable 
amount  to  occur.  One  of  the  nitrogen  atoms  transfers  suffi¬ 
cient  energy  to  the  argon  atom.  The  latter  moves  away  from 
the  bound  nitrogen  molecule  with  increased  kinetic  energy, 
as  in  Figure  10(a). 

(ii)  Trajectories  Involving  Complexes 

As  shorm  in  Figure  10(c),  one  nitrogen  atom  is  initially 
bound  to  the  argon  atom  bv  the  weak  attractive  forces,  form¬ 
ing  an  NAr  complex.  Another  nitrogen  atom  approaches  the 
complex,  and  its  strong  attractive  valence  forces  acting  on 
the  first  nitrogen  atom  pulls  the  latter  away  from  the  argon 
atom,  resulting  in  a  reactive  trajectory. 
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(111)  Trajectories  with  Multiple  Crossings  of  the  Barrier  Surface 
A  high  temperature  trajectory  of  this  type  is  shown  in 
Figure  10(b)  in  which  the  argon  atom  exchanges  energy  and 
momentum  with  each  of  the  two  nitrogen  atoms,  resulting  in 
the  reactive  transition  (b  -*■  f  ■+  b). 

At  low  temperatures,  however,  the  attractive  minimum  in  the 
NAr  interaction  potential  becomes  important,  and  a  single 
nitrogen  atom  may  exchange  energy  and  momentum  to  and  from 
the  argon  atom  more  than  once  during  a  collision,  resulting 
in  a  trajectory  with  multiple  crossings  of  the  barrier 
surface. 

Figures  11(a)  and  (b)  show  the  time  evolution  of  some  of  the  tra¬ 
jectories  examined  above.  The  upper  half  of  these  figures  plots  the 

energy  of  the  molecule  ( (H  -  B)/kT]  versus  the  time  (v  t,  where  v  is 

e  e 

the  ground  state  vibrational  frequency  of  the  molecule),  and  the  lower 
half  shows  the  quantities  and  r^  as  functions  of  time.  It  can  be 
seen  that  the  molecular  vibrations  are  highly  anharmonic,  and  the 
energy  exchanges  are  typically  of  order  kT  or  less  (see  also  Figure  12) 
For  a  collision  to  be  reactive  (consider  recombination),  the 
molecule  must  start  in  a  free  state  and  end  in  a  bound  state.  The 
corresponding  energy  changes  must  be  such  that  H  -  B  <  0,  where  H  is 
the  Hamiltonian  of  the  molecule,  and  B  is  the  height  of  the  rotational 
barrier.  Figure  12  shows  a  distribution  plot  of  AH  and  AB  for  the 
reactive  trajectories  of  the  case  N2  +  Ar  at  11500  °K  on  the  surface 
e  =  0.  In  Reference  8,  it  was  shown  that  for  HCJl  +  Ar,  recombination 
(the  change  from  H-B>0toH-B<0)  is  achieved  by  two  equally 


\ 
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lmportant  mechanisms,  i.e,,  reducing  H  and/or  increasing  B.  For 
^  +  Ar,  on  the  other  hand  (see  Figure  12),  most  recombinations  are 
achieved  by  reducing  H.  This  furnishes  further  evidence  that  rota¬ 
tional  energy  transfer  in  highly  asymmetric  molecules  plavs  a  rela¬ 
tively  important  part  in  chemical  reactions. 

The  distributions  of  the  energy  of  the  recombined  nitrogen 
molecules  are  shown  in  Figure  13  for  several  temperatures,  where  Fw 
is  the  total  energy  of  the  molecule,  is  the  height  of  the  rota¬ 
tional  barrier  (see  Reference  2),  and  (E  -  E)^  Is  a  measure  of  the 
vibrational  energy  of  the  molecule.  Figure  14  shows  plots  of  the 
distribution  of  trajectories  with  respect  to  the  Impact  parameter  (b^) 
and  the  kinetic  energy  (EL /kT) ,  respectively,  of  the  argon  atom  rela¬ 
tive  to  the  center  of  mass  of  the  nitrogen  molecule.  The  reactive 
trajectories  are  indicated  by  the  shaded  area  in  these  figures. 

For  trajectories  rrhich  involve  bound  N*Ar,  complexes,  we  define 
the  impact  parameter  (b£)  and  the  kinetic  energy  (Ec/kT)  of  the  free 
nitrogen  atom  relative  to  the  center  of  mass  of  the  N*Ar  complex.  The 
distributions  with  respect  to  these  two  parameters  are  plotted  in 
Figure  15  (for  T  ■  100  °K) .  Again,  the  reactive  (recombination) 
trajectories  are  indicated  by  the  shaded  area. 

It  is  interesting  to  note  the  characteristics  of  the  impact 
parameter  distributions.  The  number  of  trajectories  (N)  approaches 
zero  for  (b^,  b^)  ■*  0.  This  is  a  geometrical  effect.  As  (b^,  b^) 
Increases,  N  rises  to  a  maximum  and  then  dies  off  gradually  when  the 
interaction  between  the  argon  atom  and  the  nitrogen  molecule  becomes 


weaker  and  weaker 
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The  results  of  the  Monte  Carlo  trajectory  calculations  presented 
above  Indicate  that  the  equilibrium  transition  rate  for  systems  of 
diatomic  molecules  highly  diluted  in  an  inert  gas  heat  bath  can  be 
represented  satisfactorily  by  a  simple,  separable  function  which  con¬ 
tains  three  characteristic  parameters.  It  appears  that  these  parameters 
can  be  correlated  and  predicted  hv  the  formulas  obtained  in  this  study 
although  better  statistics  and  additional  data  for  other  systems  are 
probably  needed  for  certain  refinements  in  these  formulas.  Also,  in 
order  to  understand  more  about  the  effects  of  the  "complex"  mechanism 
to  the  overall  recombination  rate  coefficient  at  low  temneratures ,  a 
detailed  analysis  of  the  corresponding  master  equations  (wit^  the 
appropriate  transition  rates)  must  he  carried  out.  Finally,  it  is  of 
great  interest  to  wodify  the  present  Monte  Carlo  techniques  for  use 
with  three-body  interaction  potentials  of  a  more  general  nature:  in 
particular,  for  systems  that  involve  strongly  reacting  third  bodies, 
e.g.,  H  +  H  +  K. 
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TABLE  I 


Morse  Potential  Parameters 


System 

D(eV) 

3  (A  L) 

0 

r  (A) 
e 

N2 

9.90 

2.66 

1.10 

N2(A  V) 

3.68 

2.70 

1.29 

Mr 

0.017 

1.49 

3.18 

System  (kT/D) 


n2(x) 


(0.1) 


n2(x) 


1000 

(.0087) 


n2(x) 


200 

(.0017) 


H2(X) 


100 

(.0009) 


e 

Wi 

Do 

D2 

°4 

Do 

.210^ 

°2 

(^2  !/2 
D4 

0 

.366 

1.56 

.342 

2.14  4,26 

2.14 

1.39 

-1 

.373 

.876 

.391 

2.30 

2,34 

1.91 

1.43 

-2 

.344 

.992 

.381 

2.74 

2.83 

1.81 

1.30 

-3 

.222 

.744 

.221 

.824 

3.35 

2.01 

1.80 

mean 

3.2i 

1.97 

1.48 

9 

.501 

1.52 

.927 

17.4 

3.04 

1.08 

.80 

-1 

.395 

.600 

.988 

47.5 

1.52 

.80 

.50 

'  -2 

.280 

.479 

.587 

13.9 

1.71 

.955 

.71 

-3 

.309 

.456 

.540 

12.4 

1,48 

1.14 

.72 

mean 

ijST" 

__  - 

768 

3 

.332 

.481 

.817 

47.1 

1.45 

.813 

.46 

-1 

.229 

.234 

.732 

102. 

1.02 

.526 

.29 

-2 

.191 

.353 

.471 

14.6 

1.85 

.811 

.62 

-3 

.147 

.144 

.338 

.6.24 

.98 

.870 

.81 

mean 

1.32 

.78 

“754 

0 

.183 

.205 

.454 

11.6 

1.12 

.807 

.69 

-1 

.128 

.067 

.503 

40.2 

.524 

.509 

.39 

-2 

.102 

.053 

.438 

26.8 

.520 

.466 

.44 

-3 

.088 

.067 

.281 

16.2 

.762 

.627 

.46 

mean 

.73 

.60 

.50 
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TABLE  II 


Summary  of  Energy  Transfer  Momenta 
(Continued) 


* 

System 

T  °K 
(kT/D) 

r 

hi 

n 

D2 

»0 

mr 

2lrll 

(-J.)1/2 

D4 

Ml 

1  0 

^2 

0 

.473 

2.04 

.784 

17.5 

4.31 

1.21 

.73 

n2(a) 

1000 

-l 

.381 

.682 

.658 

9.40 

1.79 

1.16 

.92 

(.0204) 

-2 

.291 

.848 

.553 

10.4 

2.91 

1.05 

.80 

.  -3 

.233 

.606 

.387 

3.21 

2.60 

1.20 

1.20 

mean 

T790 

1.16 

.91 

0 

.360 

.782 

.792 

32.4 

2.17 

0.91 

0.54 

«2(a) 

20  0 

-1 

.291 

.314 

.886 

96.2 

1.08 

0.66 

0.33 

(.0047) 

-2 

.201 

.256 

.548 

19.9 

1.27 

0.73 

0.57 

-3 

.173 

.162 

.415 

17.9 

0.94 

0.83 

0.53 

mean 

06 

0.7$ 

0.49~ 

0 

.217 

.223 

.694 

38.8 

1.03 

.626 

.464 

k2(a) 

100 

(.0023) 

-1 

-2 

.158 

.186 

.092 

.244 

.504 

.524 

18.7 

19.0 

.582 

1.31 

.628 

.748 

.569 

.576 

-3 

.164 

.099 

.580 

34.8 

.604 

.566 

.448 

^  T  , 

mean 

"  7sW 

T42" 

”3iT 

* 

The  third  body  Is  argon. 
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TABLE  XII 


Summary  of  Characteristic  Parameters 
for  the  Transition  Kernel 


* 

System 

T  °K 

kT/D 

C 

b 

2 1 D 1 1/d2 

a 

l-(a/6)2 

11500 

.10 

.42 

.77 

1.97 

2.04 

.54 

N,(X) 

1000 

.009 

.44 

.77 

0.99 

1.12 

.76 

200 

.002 

.32 

.77 

0.78 

0.94 

.83 

100 

.001 

.25 

.77 

0.60 

0.79 

.88 

1000 

.023 

.48 

.77 

1.16 

1.27 

.72 

n2(a) 

200 

.005 

.36 

.77 

0.78 

0.94 

00 

• 

100 

.002 

.24 

.77 

0.64 

■Sir:.-: 

0.82 

.87 

* 

The  third  body  Is  argon. 
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TABLE  IV 

Distribution  of  Trajectories  with  Respect  to  Class  of 
Reaction  and  Number  of  Traverses  of  the  Energy  Surface  e. 


(a)  Ar  +  N2  at  11500  °K  (kT/D  -  0.1),  e  -  0. 


Class 

N 

Total 

n  j 

2 

3 

(f|n|b) 

103 

12 

3 

118 

(f  |  n  |  f ) 

68 

9 

1 

78 

(b  |  n  |  b) 

87 

3 

0 

90 

(b  |  n  |  f ) 

10 

3. 

1 

14 

Total 

300 

Recrossing  correction  factor 

N/Nq  -  (103  +  12/2  +  3/3)  «  0. 

37 

(b)  Ar  +  N2 

at  100  °K 

(kT/D  - 

0.0009), 

E  -  0. 

Class  n 

1  2 

N  • 

3  4” 

>  5 

Total 

(f|n|b) 

61  29 

9  4 

1 

104 

(f|n|f) 

34  9 

3  2 

1 

49 

(b  |  n  |  b  ) 

82  27 

6  0 

1 

116 

(bjnjb) 

20  8 

1  1 

0 

30 

Total 

299 

Recross ing  correction  factor 

N/Nq  -  (61  +  29/2  +  9/3  +  4/4  +  l/5)/299  -  0.27 
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TABLE  IV 

Distribution  of  Trajectories  with  Respect  to  Class  of 

Reaction  and  Number  of  Traverses  of  the  Energy  Surface  e. 

(Continued) 


(c)  Ar  +  N2  at  100  °K  (kT/D  -  0.0009),  e  -  0, 
with  non-complex  trajectories  only. 


Class 

N 

Total 

n  1 

2 

s 

“4"-' 

■"5 

(f |n|b) 

30 

15 

1 

1 

0 

47 

(f  |  n  j  f ) 

19 

3 

2 

0 

24 

(b  |  n  |  b ) 

79 

24 

4 

0 

1 

108 

(b  1  n  |  f ) 

19 

8 

1 

1 

a 

29 

Total 

208 

Recrossing  correction  factor 


N/N  -  (30  +  15/2  +  1/3  +  l/4)/208  -  0.183 
o 
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TABLE  V 


Recombination.  Rate  Coefficients  for  +  Ar 


State 

T  °K 

kT 

kB 

r 

'  N 

n“ 

o 

,  k® 
k  6 

k  ,  B 

e  kr 

«<S)(k/VEb 

D 

&(0)(k/ke)o 

"  IS. 

r 

11500 

.1 

9.92"34 

.42 

.22  .60 

^  1 

6.14"35 

X 

1000 

.009 

3. 04”  33 

.44 

.60  .77 

1.17 

7.2l"34 
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FIGURES 


Plot  of  [N(e)/N0(e)]/[1  -  exp (-  kT/D)]  versus  e  for  N2  +  Ar. 

#,  T  =  11500  °K;  ©  ,  T  *  1000  °K;  ©,  T  -  200  °K;  O  ,  T  -  100  °K. 
Solid  line  is  the  best  least-squares  fit  for  all  points.  The  slope 
of  this  line  determines  the  parameter  b. 

1/2 

Correlation  of  the  parameter  b  with  [m3/ (m^  +  m2)]  .  , 

Correlation  of  the  parameter  a  with  kT/D,  where  D  is  the  well  depth 
of  the  interaction  potential  of  the  molecule.  ▲  ,  N2(X  Ig)  +  Ar; 

A,  N2(A  V)  +  Ar;  ▼  ,  HCJl  +  Ar;  V,  +  Ar;  •,  02  +  Xe; 

O,  02  +  Ar;  ■  ,  I2  +  Ar.  Solid  line  is  the  best  least-squares 

fit  for  all  points. 

Plot  of  the  averaged  absolute  energy  transfer,  <A>,  and  the  root- 
mean-square  energy  transfer,  <A2>^2,  versus  kT/D  for  reactive 
collisions.  D  is  the  well  depth  of  the  interaction  potential  of 
the  molecule.  Remainder  of  legend  same  as  Figure  3. 

Map  of  numerical  distribution  of  the  function  Me±,  ef)  for  N2  +  Ar 
at  11500  °K.  The  contour  map,  calculated  from  Equation  (4.8) 
using  the  fitted  parameters,  is  plotted  in  the  upper  left  corner 
to  avoid  overcrowding  the  figure.  Comparison  may  be  done  bv 
reflection  through  the  diagonal 

The  function  [R(A,  e>^(0)]/lft(e>/rt(0)]  is  plotted  versus  A, 
the  absolute  value  of  energy  transfer,  for  the  case  of  N2  +  Ar 
at  11500  °K.  •  ,  e  =  0;  O,  e  -  -1;  e  -  -2;  A  ,  e  =  -3. 

The  factor  [  ^.(e)  /  <ft(0)  ]  is  included  for  normalization  purposes  so 
that  the  data  can  be  better  represented  by  a  single  line.  The 
solid  line  which  passes  through  the  bulk  of  the  data  has  a  slone 
taken  from  Table  III.  Sec  text.  Section  4.1. 
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7.  i  /”  ,  the  ratio  of  the  contribution  from  the  "complex"  mechanism 

C  PC 

to  that  from,  the  noncomplex  mechanism,  is  plotted  versus  kT/D,„j, 
where  D,„f  is  the  well  depth  of  the  third-body  interaction  potential. 


3.  Plot  of  the  effective  recombination  rate  coefficient  due  to  barrier 

penetration  versus  energy,  for  I!0  at  several  temperatures. 

9.  Comparison  between  theoretical  predictions  of  the  recombination 

rate  coefficient  and  experimental  measurements  as  a  function  of  T. 

+  ,  theoretical  predictions  of  k^(X  £u);  theoretical 

predictions  of  k  (X  1X+);  A,  Campbell  and  Thrush''1;  O  ,  Clyne  and 
r  p  ^ 

Stedtran;  -  V,  Brennen  and  Shane;  (A),  Appleton  et  al;  (B), 
Byron;35  (C),  Cary.36  Solid  and  dashed  lines,  see  text,  Section  4.2. 
10.  Typical  trajectories  for  +  Ar  collisions. 

(a)  Shows  f  -»•  b  transition  at  11300  °K. 

(b)  Shows  b  -*  f  -►  b  transitions  at  11500  °K.  O  »  position  when  b  -*•  f 
transition  occurred;  0,  position  when  f  b  transition  occurred. 

(c)  Shows  f  -*■  b  transition  at  100  °K.  In  the  upper  part,  the  Ar 
atom  is  at  the  origin;  whereas,  in  the  lower  part,  the  origin 
corresponds  to  the  center  of  mass  of  N0.  The  motion  of  the 
l!*Ar  bound  complex  is  clearly  illustrated. 

The  small  dots  on  the  trajectories  indicate  time  intervals  of  Vgt  =  1 
in  (a)  and  (b)  and  Vgt  =  10  in  (c)  where  Vg  is  the  ground  state 
vibrational  frequency  of  the  molecule. 
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11.  Time  evolution  of  typical  trajectories.  Upper  half  of  the  figure 

shows  energy  of  the  molecule  [(H  -  B)/kT]  versus  time  (V^t). 

Lower  half  shows  r,  „  and  r»  versus  V  t.  r  is  the  equilibrium 
12  3  e  e 

separation,  and  Vg  is  the  ground  state  vibrational  frequency  of 
the  molecule. 

(a)  N2  +  Ar  at  11500  °K  , 

(b)  N2  +  Ar  at  100  °K  . 

12.  Distribution  plot  of  reactive  trajectories  in  the  [AH,  AB]  plane 
for  N2  +  Ar  at  11500  °K. 

13.  Energy  distributions  of  nitrogen  molecules  formed  by  recombination 
in  argon. 

14.  Distributions  of  trajectories  (a)  with  respect  to  the  impact 
parameter;  (b)  with  respect  to  the  kinetic  energy  of  the  third 
body  (see  text.  Section  V)  for  N2  +  Ar  at  several  temperatures. 
Reactive  trajectories  are  indicated  by  the  shaded  areas. 

15.  Distributions  of  trajectories  with  respect  to  the  impact  parameter 
and  the  relative  kinetic  energy  defined  for  the  complex  mechanism 
(see  text,  Section  V),  for  N-Ar  +  N  at  100  °K.  Reactive  (recombina¬ 
tion)  trajectories  are  indicated  by  the  shaded  areas. 
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